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Abstract 

We construct a sequence of functions that uniformly converge (on compact sets) to the 
price of Asian option, which is written on a stock whose dynamics follows a jump diffusion, 
exponentially fast. We show that each of the element in this sequence is the unique classical 
solutions of a parabolic partial differential equation (not an integro-differential equation). 
As a result we obtain a fast numerical approximation scheme whose accuracy versus speed 
characteristics can be controlled. We analyze the performance of our numerical algorithm 
on several examples. 

Key Words. Pricing Asian Options, Jump diffusions, an Iterative Numerical Scheme, 
Classical Solutions of Integro-differential equations. 

1 Introduction 

We develop an efficient numerical algorithm to price Asian options, which are derivatives whose 
pay-ofF depends on the average of the stock price, for jump diffusions. The jump diffusion 
models are heavily used in the option pricing context since they can capture the excess kur- 
tosis of the stock p rice returns and along with the skew in the implied volatility surf ace (see 



tosis oi tne stocK price returns and along witn tne sRew m tne impued voiatiiity surt ace (see 
Cont and Tankovl (|2003l l). Two well-known examples of these models are i) the mo del oflMertonl 
m which the jump sizes are log-normally distributed, and ii) the model of iKoul (j2002l ) 
in which the logarithm of jump sizes have the so called double exponential distribution. Here 
we consider a large class of jump diffusion models including these two. 

The pricing of Asian options is complicated because it involves solving a partial differential 
equation (PDE) with two space dimensions, one va riable a ccoun ting for the averag e stock 
price, the other for the stock price itself. However, Vecer ( 2001 ) and Vecef and Xu ( 20041 ) 
were able to reduce the dimension of the problem by using a change of meas ure argumen t 
(also see Section 2.1). When the stock price is a geometric Brownian motion IVeceii ( 200ll ) 
showed that the price of the Asian option at time t = 0, which we will denote by — > V{So), 
satisfies V{Sq) = Sq ■ v{z = z*,t = 0) for a suitable constant z* , in which the function v 
solves a one dimensional parabolic PDE. When the s tock price is a jump d iffusion, then under 



the assumptions that vt, Vz and Vzz are continuous Vecef and Xu ( 20041 ) (Theorem 3.3 and 



Corollary 3.4) showed that the function v solves an integro-partial differential equation using 
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Ito's lemma. However, a priori it is not clear that these assumptions are satisfied. In this 
paper, we show that for the jump diffusion models these assumptions are indeed satisfied (see 
Theorem 12. ip . i.e. we d irectly show th a .t v is the unique classical solution of the integro-partial 
differential equation in IVecef and hmi \. We do this by showing that v is the limit of a 
sequence of functions constructed by iterating a suitable functional operator, which we will 
denote by J, that takes functions with certain regularity properties into the unique classical 
solutions of parabolic differential equations and gives them more regularity and finally showing 
that V is the fixed point of the functional operator and that it satisf ies the certain regularity 
properties. This proof technique is similar to that of iBavraktaii (j2007l ). in which the regularity 
of the American put option prices are analyzed. In the current paper, some major technical 
difficulties arise because the pay-off functions we consi der are not b ound ed and also because the 
sequence of functions constructed is not monotonous ([Bavraktarl (|2007l ) was able to construct 
a monotonous sequence because of the early exercise feature of the American options). 



The iterative construction of the sequence of functions which converge to the Asian option 
price naturally leads to an efficient numerical method for computing the price of Asian options. 
We prove that the constructed sequence of functions converges to the function v uniformly 
(on compact sets) and exponentially fast. Therefore, after a few iterations one can obtain the 
function v to the desired level of accuracy, i.e. the accuracy versus speed characteristics of the 
numerical method we propose can be controlled. On the other hand, since each element of 
the approximating sequence solves a parabolic PDE (not an integro-differential equation), we 
can use one of the classical finite difference schemes to determine it. We propose a numerical 
scheme in Section [3] and analyze the performance it in the same section. 



So far nunierical methods f or the pricing Asian options were proposed for diffusion models: 
Vecefl (|2nnih . [Rogers and Shil (Il99,^ ^. Izhansl \2mi \ de veloped PDE me thods, iGeman and Yon 



(j 19931 ) developed a single Laplace inver sion method, iLinetskv I (jiooi) developed a spectral 
expansion appr oach, Cai and Koul (^ 20071 ) d eveloped a d ouble Laplace inversion method. On 
the other hand Rogers and Sh i (19 9^) and ThomDsonI ( 19981 ) obtained tight bounds for the 
prices. We should mention that ICai and Kou (2007) also considered pricing Asian option for 
the double exponential jump diffusion model of iKoul (120021 ). 



The rest o f the p aper is organized as follows: In Section 2.1, we summarize the findings of 



Vecef and Xul (|2004l ) in the context of jump-diffusion models. In Section 2.2, we present our 



main theoretical results: Theorem 12.11 and Corollary 12. li We propose a numerical algorithm 
in Section 3.1 and analyze the convergence properties of our algorithm in Section 3.2 (see 
Propositions 13. Il and [3T2|) . Then we perform numerical experiments for two particular jump 
diffusion models and analyze the performance of our numerical algorithm in Section 3.3. Section 
4 is devoted to development of the proof of Theorem l2.ll In Sections 4.1 and 4.2 we develop the 
properties of the functional operator J and the properties of the sequence obtained by iterating 
J, respectively. These results are used to prove Theorem 12.11 in Section 4.3. A brief summary 
of our proof technique can be found in Section 2.2. right after the definition of the operator J 
in (f2TT]) . 
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2 A Sequential Approximation to Price of an Asian Option 



2.1 Dimension Reduction 



Let (ri, T ^ P) be a complete probability space hosting a Wiener process {Bt;t > 0} and a Poisson 
random measm'e N, whose mean measure is \v[dy)dt, independent of the Wiener process. Let 
i^t)t>o denote the natural filtration of B and N. In this filtered probability space, let us define 
a Markov process S = {St; t > 0} via its dynamics as 



dSt 



(r - i^)St- dt + a St. dBt + St- ! {y-l 



)N{dt,dy), 



(2.1) 



in which r is the risk free rate, fi = A(^ — 1) with assumption ^ = J^^ yv{dy) < oo. The 
process S is the price of a traded stock, and under the measure P, the discounted stock price 
[e~^^ St) is a martingale. In this framework the stock price jumps from at time t the stock 
price moves from St- to St-Y, in which y's distribution is given by u. 1" is a positive random 
variable and note that when Y < 1 then the stock pric e S .jumps dow n when Y > 1 the stock 
price jumps up. In Merton's jump diffusion mo del (s ee MertonI ( 19761 )) Y = exp(X) where X 



is a Gaussian random variable. In Kou 's model (|KouI ^200i )) Y = exp{X) in which X has the 
double exponential distribution. 



To reduce the dimension of the Asian option pricing problem. IVecef and Xul (|2004l ) introduce 
a new measure Q by 



dl 



e-rt^^ tG[0,T]. 
bo 



(2.2) 



Here, T is the maturity of the Asian option. IVecef and Xul (|2004 ) also introduce a numeraire 
process 



J A 



Xt 



t tG[o,r], 



(2.3) 



where X = {Xt] t G [0, T]} is a self- financing portfolio, which replicates the pay-off of the Asian 
option, whose dynamics are given by 



dXt = qtdSt + r{Xt- - qtSt-) dt, Xq 



X, 



in which qt defined as 



t G [o,r]. 



is the number of shares invested in the stock, and 

X = qoSo - e 



-rT 



(2.4) 
(2.5) 

(2.6) 



Vecef and Xul (j2004l ) showed that the price of the continuous averaging Asian option with 



floating strike Ki and fixed strike K2 defined by 



-rT 



Stdt-KiST-K2 







(2.7) 
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can also be represented as 

V{So) = So-El^A{C-{Z^-K,))\ (2.8) 

in which G {—1, 1} indicates whether the option is a put or a call. Throughout this paper, 
the short hand notation E^^ represents the conditional expectation under Q, given the process 
at time t is z. Under the measure Q, the Poisson random measure will have mean measure 
XD{dy)dt with new jump measure = yv[dy). 



2.2 Main Theoretical Results 

In this section we will show that 

V{So) = S^-v{Zi,Q), (2.9) 

for some v that is the classical solution, i.e. v G C^'^, of an integro-partial differential equation 
(see Theorem 12.11 and Corollary 12. ip . We will also show that v is the limit of a sequence of 
functions constructed by iterating a functional operator, which is defined in (j2.1ip . We will 
show that each of the functions in this sequence are classical solutions of partial differential 
equations (not integro-differential equations) and that they converge to v locally uniformly 
and exponentially fast (see Theorem 12. ip . The analytical properties of the functional operator 
(listed in the lemmas of Section I4.ip used in the construction of the approximating sequence 
play important roles in proving our main mathematical result. We will summarize the role of 
the functional operator below after we introduce it. 
Let us introduce the following sequence of functions 

vo{z,t) = {C- {z- Ki))-^ , Vn+i{z,t) = Jvn{z,t) n > 0, for all (z,t) G M X [0,r], (2.10) 

in which the functional operator J is defined, through its action on a test function / : Mx [0, T] — > 
M-i-, as follows: 



j/(z,t) = e; 



t,z 



|e-A«(T-t) . - Ki))+ + e-^«(^-*) A • P/(Z„ .) d.j , (2.11) 



in which Z = {Zt;t >0} has the dynamics 

dZt = -fi {qt - Zt) dt + a {qt - Zt) dWt 

and 

PfiZt,t) = [ f(zt + {qt - Zt)^, t) yv{dy) 

) ' ^ (2.12) 

= / fl^+qt'^Ay.idy). 

Jr+ \y y J 

We will show that the sequence of functions defined in (|2.1Up by iterating J are classical 
solutions of PDEs thanks to the following analytical properties of the operator J (which are 
developed in Section 14. ip : 1) J maps functions that are Lipschitz continuous with respect to 
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the z- variable (uniformly in the variable) and Holder continuous with respect to the t- variable 
into classical solutions of PDEs (see Proposition I4.ip . 2) J preserves Lipschitz continuity with 
respect to the z-variable (see Lemma [4.ip . 3) J transforms Lipschitz continuous functions, with 
respect to the z-variable, that satisfy a linear (in the z-variable) growth condition (uniformly 
in the t- variable) into Holder continuous functions of the t- variable (see Lemma 14. 3p . 4) J 
preserves the linear growth condition in the z— variable (see Lemma 14.21 and Remark 14. ip . The 
analytical properties of J can be summarized as "J maps nice functions (set of functions with 
a few regularity properties), to nicer functions (set of functions that are the classical solutions 
of partial differential equations, and have the same regularity properties as before). 

It is a priori not clear that the sequence of functions defined in ()2.10p has a limit. Using 
the properties of the operator J we show that this sequence is Cauchy (see Lemma 14. 6p and 
therefore has a limit (in fact the sequence converges locally uniformly and exponentially fast). 
We show that, the limit of this sequence, which we denote by Voo^ is a classical solution of an 
integro-PDE using 1) the fact that it is a fixed point of the operator J (see Lemma l4.7p . 2) the 
facts that it is Lipschitz continuous with respect to the z- variable (uniformly in the variable) 
(see Lemma l4.8p and Holder continuous with respect to the t-variable (see Corollarv l4.3p . 

Finally, using a verification argument we will show that the limit Voo is indeed the function 
that satisfies (j2.9p (see Corollarv l2.ip . 

The main theoretical results that are summarized above will be stated in the next theorem 
and its corollary. The proof of Theorem 12.11 is given in Section 14.31 which uses the results of 
Sections 14.11 and 14.21 as summarized above. 

Theorem 2.1. (i) The sequence of functions defined in i2.10\) has a pointwise limit. Let us 
denote this limit by Voo{z,t). 

(ii)For any compact domain P C M, Vn{z,t) converges uniformly to Voo{z,t) for {z,t) G 
V X [0,T]. Moreover, 




(2.13) 




(2.14) 
(2.15) 



At) = -Kqt - z)- + -a\qt - zf^. 




(2.16) 



(iv) The function is the unique classical solution, i.e. Voo G C^'^, of 



A{t)Voo{z, t) - X^Vooiz, t) + X - {PVoo){z, t) + -^V^{Z, t) = 

vUz,T) = {C ■ {z - K^))+ . 



(2.17) 
(2.18) 



Proof. See Section 14.31 



□ 
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The iterative procedure in (j2.14p simply collapses to the Vecer's PDE (see Vecer ( 200ll )) 
when A = 0, i.e. when the underlying asset is a geometric Brownian motion. That is, the 
iteration in (j2.14p is designed for the models in which the asset price jumps. 



Corollary 2.1. Let V{So) be as in {2.1), i.e. V{So) is the value of the Asian option for jump 
diffusion S whose dynamics is given in \2.1\) . Then we have 



V{So) = So ■ Vooiz,0), 



(2.19) 



in which foo(-,-) is the unique solution of the integro-partial differential equation ^2.11^ with 
terminal condition h2.18\) . and 



Proof. Let us define 



(2.20) 



(2.21) 



where Z"^ , defined in (12. 3p . has the initial value = z. It follows from (12.170 and the Ito's 
lemma that Mt is a Q-martigale, i.e. Mt = E'Q{Mt|J^J. As a result 

v^{z,Q) = Mo = eJJMt} = Egj^;oo(^^,r)} = Egj(C • {Z^ - Ki))+} = (2.22) 



The last identity follows from the representation ()2.8p . 



□ 



3 Computing the Prices of Asian Options Numerically 



It follows from Theorem 12.11 that the the sequence of functions (wn)n>o converges uniformly 
and exponentially fast to Voo on any compact domain. Therefore a few iterations of (j2.14p 
and (j2.15p . starting from vq will produce an accurate approximation to Voo- To perform the 
iterations we will make use of the finite difference methods for PDEs. Since each Vn+i{-, ■) is the 
classical solution of a partial differential equation (not an i ntegro-partial d i fferent ial equation) 
we can use Crank-Nicolson discretization (see page 155 of Wilmott et al. ( 19951 )) along with 
the SOR algorithm (see e.g. page 150 of Wilmott et al.l ( 1995 )) to solve the sparse system of 
linear equations. In this iteration we will need to compute the integral Pvn, and we do this by 
the trapezoidal rule. 

We will describe this numerical method more precisely in Section 13.11 and investigate the 
convergence properties in Section [3.21 In Section [3.31 we determine the performance (the s peed 
and a ccura cy characterist ics) of our numerical method for the jump diffusion models of IKou 
(|2002l ) and lMertonI (|l97fil ). i\ this section we will take the Monte-Carlo simulation results as a 
benchmark. 



3.1 A numerical algorithm 

Let us discretize (j2.17p using the Crank-Nicolson method. For fixed At, Az, 2;jnax and Zmin, let 
MAt = T and KA.z = Zmax — -Zmin- Let us denote Zk = Zmin + ki^z = 0, 1, • • • ,K. By v we 
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will denote the solution of the difference equation 



(1 + Pk,m)vik, m) - pt^„,v{k + 1, m) - Pi, .^v{k - 1, m) 

= Pk,m+i^{k + 1, m + 1) + Pk^^_^iv{k - 1, m + 1) + (1 - pl^^+i)v{k, m+l) 



Pv] {k,m + I) + [Pv] (k, m) 



(3.1) 



for m = M — 1, M — 2, •, 0, k = 0, 1, • • • , i^T, satisfying the terminal condition v{k, m) = {C ' 
{kAz — Ki))~^ and appropriate boundary conditions. The coefficients p^^, P^m ^^'^ Pkm 
given by 



PL 



Pk, 



a 



a 



2 ( q niAt - kAz 
Az 

2 ( Q mAt - kAz 

Az 



_ Qm At - kAz 
Az 

QmAt - kAz 

Az 



At, 



At, 



(3.2) 



Pk,m=Ptm+Pk,m + 2^^^^- 



In (j3.ip . P is a linear operator which is the discrete version of the operator P in (j2.12p . 
Letting x = log y, we can write Pf as 



f 



z e - 



1 



t e''F{dx) 



(3.3) 



in which F{dx) is the distribution of a random variable X, such that the distribution of 
is given by the jump measure v. We approximate the integral in (j3.3p using trapezoidal rule. 
Discretizing a sufficiently large interval [xmin-,Xmax\ into L subintervals, we obtain the grid 
Xmin = xq < xi < • • • < XL = Xmax- This grid may not be equally spaced. One can choose the 
grid to be finer where density of the distribution F is large. The left hand side of ()3.3p will be 
evaluated on grid point {zk,mAt). However for some m, k and £, z/^/e^^ + qmAt{e^' — l)/e^^, 
as the first variable of / in the integrand, may not land on z^/ for some k' . Consequently, we 
will determine the value of the integrand in (j3.3p by linear interpolation. If 

^k e^* - 1 , , , 

Zk' < + QmAt z— < zyj^i, for some k 



then 



f (— + qmAt—^,mAt 



(1 - w)f{zk',mAt) + wf{zk'+i,mAt) + 0((Az) 



in which w is the interpolation weight. On the other hand, if z^je^^ + QmAtis^' — l)/e^^ is 
out side the interval [zmin, z:max], the value of the function is determined by the boundary 
conditions. Now the integral in (j3.3p can be approximated as 
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f(^ + gt'^,t]e^Fidx) 



1 r / 



£=0 



+ QmAt—^^ — ,mAt e^^+i5r(x£+i) 



{xe+i-xe)+OiiAxf), 



(3.4) 



where Ax = max^ — a;^), and g is the density of F. 

Note that numerically solving the system of equations in (13. ip is quite difficult due to the 
contributions from the integral terms (i.e. the Pv{k,m) term). Discretizing ()2.14p recursively 
(using the Crank-Nicolson discretization) we obtain sequence Vnik, m), = 0, 1, • • • , by setting 
VQ{k,m) = (^ • (kAz — i^i))"*" through the recursive relationship 

(1 + Pk,m)vn+i{k,m) - p+^^Vn+i{k + l,m) - p-^^^^Vn+iik - l,m) 

= Ptm+l'^n+l{k + l,m+l) +p';^^^^Vn+l{k-l,m + l) + (1 -p2,m+l)^'n+l(^,"l + 1) 



+ -XAt 



PVr 



m + 1) + [Pvn] {k,m 



(3.5) 



For each n the terminal condition Vn+iik,m) = ((" • {kAz — i^i))"*" and appropriate boundary 
conditions are s atisfied. We wi l l solv e the sparse linear system of equations using the SOR 
method (see e.g. Wilmott et al. ( 19951 )). 



3.2 Convergence of the Numerical Algorithm 

In what follows we will first show that as n ^ oo, Vn converges to v. Next, we will argue that 
as the mesh sizes At and Az go to zero v converges to Voo- For the sake of simplicity of the 
presentation, in what follows we will assume that {P\){k,m) < j^^yv{dy) = ^. Otherwise 
the order of error of the discretization of the integral will have to be sufficiently small for the 
following statement to be true. 

Proposition 3.1. For sufficiently small At and Az, Vn converges to v uniformly and at an 
exponential rate. 

Proof. Let us define 

en{k,m) = v{k,m) — Vn{k,m). (3-6) 

Since P is a linear operator e„ will satisfy (j3.5p when we replace Vn by and Vn+i by e^+i. 
Now let us define 

E"^ = max |e„,(A:, m)|, £'n = max |e„,(/c, m)|. (3-7) 

k m,k 

Choosing At and Az small enough we can guarantee that p^^, P^m ^ ~ Pkm positive 
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for all {k,m). As a result it follows from the difference equation e„ satisfies that 



(1 + pU K+i{k, m)| < (p+„^ + + {plm+i + Pk,m+i + (1 - Plm+i))K-^i + AAtC^„ 

(3.8) 



in which we used the assumption that {Pl){k,m) < S^. It follows from (j3.8p that 

(1 +pO,„J |e„+i(/c,m)| - (pO,^^ - ^AeAt) i^™ ^ < (l - lA^At) + AAt^i^^. (3.9) 

Let A;* be such that |e„+i(A;*, m)| = El^;^-^^. Since the right-hand-side of (|3.9p does not depend 
on fc, we can take k = k* on the left-hand-side to write 

£^™+l<eKVY + (l-^)^n, (3.10) 

in which 

1 + 1/2 • A^At ^ ' ^' ^ ^ 

as a result of the assumption p^^, P^^ ^ ~ Pkm positive for all {k, m). It follows from 
([3111]) that 

< e^'-'^E^'^^ + (1 _ + + . . . + 0M~m-l^En. (3.12) 

Because of the terminal condition of Vn, we have -E^'^i = 0. In addition, (j3.12p is satisfied for 
all m we get that 

En+l < (1 - 0^')En. (3.13) 

As a result 

< (1 - 0^^)'^Eo -^0 as n ^ oo. (3.14) 

□ 

Remark 3.1. As M ^ oo 

which shows that the convergence rate in ^3.14^ agrees with the convergence rate in \4.6l ). 
Proposition 3.2. 

\voo{zk-,rnAt) — v{k,m)\ ^ Q (3.16) 

as Az, At, Ax 0. 

Proof. Using the triangle inequality let us write 

\voo{zk,mM) - v{k,m)\ < \voDizk,mAt) - Vnizk,mAt)\ + \vn{zk,mAt) - Vn{k,m)\ 

+ \vn{k, m) — v{k, m)\ 

<c(l- e-^''(^-™^*))" + n • 0{{Atf + {Azf + (Axf) + C(l - 0^^)", 

(3.17) 
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for some positive constants C and C . The first and the third terms in the right-hand-side of the 
second inequahty are due to (I4.61|) and ()3.13p . The second term arises since the order of error 
from discretizing a PDE using a Crank-Nicolson scheme is 0((Az)^ -|- (At)^), the interpolation 
and the discretization error from the numerical integration are of order (Az)^ and (Ax)^ and 
that the total error made at each step propagates at most linearly in n when we sequentially 
discretize the PDEs in ([2T4]) . 

Letting Ai, Az ^ in (j3.17p we obtain that 

^^hm ^\v^{zk,mM)-v{k,m)\<c(l-e-^'''^Y + c{l-e-^^'^Y , (3.18) 

in which we used (j3.15p . Since n is arbitrary the result follows. □ 

Remark 3.2. The proof of Proposition \3.^ leads itself to an order of convergence. In \3.11^ 
we can choose n = 0(log(l/[(At)^ + (Az)^ -|- (Ax)^])), which would guarantee that the first 
and the third terms on the right-hand- side are of order 0((Ai)^ -|- (Az)^ -|- (Arc)^), and the 
right-hand- side becomes of order [(At)^ + {/\zf + (Ax)^] log(l/[(At)2 + (Az)^ + (Ax)^]). Note 
that this order of convergence is better than that of 0((At)^~'>' + (Az)^"'*' -|- (Ax)^"'*') for any 
7 > 0. 



3.3 Numerical results for Kou's and Merton's models 



There are t wo we l l-know n examples of jump diffusion in the li t eratu re, the double exponential 
model as in IkouI \2m± and the normal model as in iMertonI (Il976l ). In this section, we will 
demonstrate our algorithm in Section 13.11 in pri cing Asian optioi is for th e se tw o models. We 
will introduce the jump distributions chosen by iKou (|2002l l and iMertonI (|l976l ^ next. Let X 
be a random variable whose probability distribution function is equal to a given distribution F 
and let the jump measure v be equal to the distribution of the random variable e^. In Kou's 
model, F is the double exponential distribution whose density is 



F{dx) = (p??ie-^^n|,>o} + (1 -p)r?2e^^n{,<o}) dx. 



In Merton's model, F is the normal distribution whose density is 

1 f-{x-iif 



F{dx) 



V27r(5-2 



exp 



^2 



dx. 



(3.19) 



(3.20) 



The price of the Asian option with floating strike Ki and fixed strike i^2; whose pay-off 
function is given by (j2.7p . can be calculated in terms of f 00(2,0) as a result of Corollary 12.11 
On the other hand, in Proposition 13.21 we have shown that Voo can be approximated by its 
discrete version v. Furthermore, Proposition 13.11 tells us that v can be approximated by the 
sequence {vn)n>o with exponentially fast convergence. Therefore, we will approximate the price 
of Asian options by iteratively solving a sequence of difference equations in (j3.5p . which are 
discretization of a sequence of parabolic partial differential equations (not integro-differential 
equations) given by (|2.14p and ()2.15p . 

In the following, we will list the numerical results for the prices of Asian options. Since 
we could not find any numerical results on European Asian options for jump diffusion models 
in the literature, we use the Put-Call parity for Asian options as a consistency check for our 
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results. For the European Asian option with floating strike Ki and fixed strike K2, the put - 
call parity gives the following identity between call and put option price 

C{So,0) - P(5o,0) = ^(1 - e-^'^)5o - K^So - e-^'^ K^. (3.21) 
ri 

This identity does not depend on dynamics of the underlying process S. Using our algorithm, we 
will calculate the call and put option price independently. Then we will compare the difference 
between our call and put price with the difference predicted by (j3.2ip . In addition, we will also 
compare our numerical results with the Monte Carlo results. 

In Tables 1 and 2, we will list the numerical results for the prices of European Asian options 
for both Kou's model and Merton's model. Run times are in seconds. All our computations are 
performed on a Pentium IV 3.0 GHz machine with C++ implementation. In Tables 3, 4 and 5, 
we will list the convergence results of the for the Asian option prices for the double exponential 
jump model. The parameters for both call and put options are the same as the 7th row in 
Table 1, i.e. r = 0.15, Sq = 100, Ki = 0, = 90, T = 1, a = 0.2, A = 1 and 771 = r/2 = 25. 

As we can see from these tables, our algorithm is stable with respect to all parameters and 
the convergence is fast. Moreover, our difference between call and put option prices are within 
±0.01 comparing to the difference predicted by put-call parity in ()3.2ip . The call option prices 
are almost within the standard error of the Monte Carlo results. 

If the dynamics of Z"^ only contains the diffusion part (i.e. A = 0), our algorithm is simply 
SOR. Using the same parameters chosen in IVecefl (l200lh . the SOR gives approx i mate option 
price with error ±1 x 10~^ dollar (we compared our results to Table 2 in IVecerl jMll)), and 
the run times are below 0.02 second. Comparing the numerical results in Tables 1 and 2, we 
see that the evaluation of the integral term by numerical integral is the time consuming part. 
This c an be speeded up by using the Fast Fourier Transform (see e.g. lAlmendral and Qosterlee 
(l2007l ) for the application of the Fast Fourier Transform in American options under the Variance 
Gamma model). 



4 Mathematical Analysis 

The purpose of this section is to provide the necessary background to prove Theorem l2.1[ First, 
in Section 14.11 we will analyze the properties of functional operator J: We will analyze how J 
increases the regularity of certain class of functions. J takes functions with certain regularity 
properties into the unique solutions of parabolic differential equations and gives them more 
regularity (see Lemmas 14. 1114.31 and Proposition 14. ip . 

Next, in Section [4.2^ we will develop the properties of the functions defined in (j2.10p in a 
sequence of lemmas and corollaries using the results developed in Section [4. 11 These properties 
will then be used to prove Theorem 12.11 in Section [4.31 



4.1 Properties of Operator J 

First, we will develop a representation of the functional operator J that is ame nable to regul arity 
analysis, which is carried out in this section. Using the notation in page 8 of Pham ( 19981 ). we 
can rewrite J as 

cT-t 



J/(z,t)=E« e-^«(^-*)(C-(4l,-Ki)) + / e-^«^A•P/(Z*.^^ + .)d. 



(4.1) 



11 



in which 



Pf{Zi'',t + s)= fi^ + q,^,^y^,t + s]yu{dy). (4.2) 

Jr+ \ y y J 

the process Z*'^ = {Z^'^; s > 0} has the dynamics 

dZ'/ = -fi{qt+s-Zl'')ds + a{qt+s-Zl'')dWs, Z^^' = z, (4.3) 

where {VFs}s>o is a Wiener process under the measure Q. It is possible to determine the 
solution to (14. 3p explicitly. For this purpose it will be convenient to work with the process 
Zs — qt+s — Zl'^. It follows from (j4.3p that the dynamics of Z are given by 



dZs = fiZsds — aZgdWs + g{t + s)ds, Zq = z = qt — z, (4-4) 

in which g{t) = ^qt- Now it is easy to obtain the solution of stochastic differential equation 
(lOll as 

Zs = zH^+ H^_^g{t + v)dv for s > 0, (4.5) 
Jo 



in which 



H^AeM{f^-l^^)s-aWs). (4.6) 



As a result we have that the solution of (j4.3p is given by 

Z'/ = zH^^ +bs, s>0, (4.7) 

in which 

bs ^ qt+s - qtH^s - r 9{t + v)dv. (4.8) 
JO 

It follows from ()4.7p that the solution of the stochastic differential equation ()4.3p is linear with 
respect to its initial value z. Inserting its solution (j4.7p back into the definition of the operator 
J in ()4.ip . we obtain 

Jf{z, t) = E^^ |e-^«(^-*) (C • [zH^T-t + bT-t -K,))+ + iy e'^^^A • PfizH^ + 6„ t + s) dsj . 

(4.9) 

In the following, we will study the regularity properties of the operator J with respect to 
both space and time. When the function / is Lipschitz continuous with respect to its first 
variable, the following lemmas show Jf is not only Lipschitz with respect to its first variable, 
but also Holder continuous with respect to the second variable. 

Lemma 4.1. For any t € [0, T], let us assume the function f satisfies 

\f{z,t)- fiz,t)\<D\z-z\, z,zGR, (4.10) 
for a positive constant D that only depends on T . Then Jf satisfies 

\Jf{z,t) - Jf{z,t)\ < E\z- z\, z,z£R, (4.11) 

with E = max{l, D}. 
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Proof. ^From the definition of operator J in (j4.9p . we have 

\Jf{z,t) - Jf{~z,t)\ < E<^|e-"«(^-*) |(C • {zH^.t + bT-t - K,))^ - (C • {zH^, + hr-t - K,))' 



T-t 



+ / dse-^^'\-\Pf{zH^^+bs,t + s)-Pf{zH^^+bs,t + s)\ 



(4.12) 



Let us obtain a bound on the right hand side of (|4.12p . First observe that 

\{C ■{zH^T-t + bT-t - Ki))+ -{C-{zH^T-t + bT-t- Ki))+\ < \z - ~z\H^_t, (4.13) 



and 



\Pf{zH^ + bs,t + s) - Pf{zH^,+b,,t + s) 



< 



f 



zH2 bs y-l , 

H ^qt+s — ,t + s ] - f 



y y 



y 



fzH^ bs 

\ y y 



t + s 



< [ D\zH^,-zH^,\u{dy) 



D\z-z\H'^s. 



yv{dy) 



(4.14) 



On the other hand, from the definition of in (j4.6p . we have that 
Inserting (|4.13p . (|4.14p and ()4.15p back into the equation (j4.12p . we have 



(4.15) 



\z — z\ e 



,(M-AO(r-t) 



fT-t 

\Jf{z,t)- Jf{z,t)\ < e-^^^^-'^z - z\E'^{H^_t} + / dse-^^' XD\z-z\E^{H^} 

Jo 

.T-t 

Jo 

< (^z) + (i_i))e-A(T-t)' 

< max{l, D} \z — z\. 



\z — z\ 



(4.16) 
(4.17) 

□ 



Remark 4.1. Let us define 



Mf ^ sup /(0,t), 

te[0,T] 

Mjf ^ sup Jf{0,t). 

te[o,T] 



(4.18) 
(4.19) 



It follows from the Lipschitz conditions ( [^. jO[ ) and ( [^. iip that both f and Jf satisfy linear 
growth conditions, if Mf and Mjf are finite, since for {z,t) G M x [0,T] 



f{z,t) < f{0,t) + D\z\, 
Jf{z,t) < Jf{{),t)+E\z\. 



(4.20) 
(4.21) 
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In the next two lemmas we will need the following moment estimates of Z*'^. 



]EQ||^M|| < C{l + \z\), 
E^{\Zl''-z\} < C{l + \z\)s^, 



(4.22) 
(4.23) 



in which < s < T and C is a constant depending on T. These estimates can be found in 
PhamI (|l998l ) (Lemma 3.1). 



Lemma 4.2. We have that 



Mjf <U + a{Mf + j] , 



(4.24) 



in which a = 1 — e ^^"^ < 1, and U, B are positive constants depending on T. 

Proof. We will estimate Mjf using the definition of the operator J in (|4.1|) . First, we have that 



C-iz. 



7tfi 



-X^T-t) 



7tfi 



in which we obtain the last inequality using the expression of Z*'^ in (jiTTl) with z = 0. First, it 
follows from (14.220 with z = that 

E'^{\bT-t\} = J} < C. (4.25) 

Letting U = C + Ki, which is a finite positive constant depending on T, we have that 

EQ |e-AC{T-t) +Ki)^< U. (4.26) 

Second, we will estimate the second term in the definition of J in (|4.ip . From the definition 
of Pf in ()4.2p . we have 



{Pf{Zl'^,t + s)} 



/ 



Z..' y-l 



+ Qt+s 



,t + s] yu{dy) 



Z. 



+ D qt+s 



y-i\ 



< E^< / I f{0,t + s) + D 

y ' y 

< e/(0, t + s)+D qt+s{^ + l) + DE'i{\ Zf \ ] 

< Cf{0,t + s)+Dqt+s{C+l) + C-D. 



yv{dy) 



(4.27) 



To obtain the first inequality we use the inequality (j4.20p , whereas the second inequality follows 
from ly — 1| <y + l. To obtain the last inequality, we use the inequality ()4.22p with 2 = 0. 
Now, using ()4.27p we obtain 



E"^ 



T-t 



e-^^'X- Pf {Zl'°,t + s) ds \ < 



T-t 



e-^«^A-[e/(0, t + s)+D qt+s{^ + I) + C ■ D] ds. 

(4.28) 
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Since < s < T — t, we have qt+s < f^- Let us define 

1 



B 



rT 



+ + c 



(4.29) 



which is a finite positive constant depending on T. Now, we have the following estimation on 
the left hand side of (|i:28]) 



T-t 



T-t 



e-^^'X- Pf{Zl'°,t + s)ds} < / e-^^'X-[U{0,t + s) + B]ds 



B 



^Prom inequalities (|4.26p and (I4.30p . we conclude that 

J/(0,t)<C/+(l-e-^«^) (^/ + |)- 



Remark 4.2. Lemma \4.^ and Remark \4-1\ indicate that 

f{z,t) < Mf + D\z\<D{l + \z\), 
Jf{z,t) < U + a(^Mf + j^+E\z\<E{l + \z\ 



(4.31) 
□ 

(4.32) 
(4.33) 



in which D = max{Mf,D} and E = max{[/ + a{Mf + B/(^),E}. We will use these linear 
growth properties to show a regularity property of the operator J with respect to time in the next 
lemma. 



Lemma 4.3. Assume the function z ^ fiz,t) satisfies 

\f{z,t)-f{z,t)\<D\z-z\, 
for z, z & W as in Lemma \4Tl\ and Mj < oo. Then 1 1— > Jf{z, t) satisfies 

\Jf{z,t)-Jf{z,s)\<F{l + \z\){s-t)K 0<t<s<T, 
in which F is a positive constant that only depends on X, ^, T and Mf. 



(4.34) 



(4.35) 



Proof. For any h G \t,T], it follows from the definition of operator J in (j4.ip and the Markov 



property of Zl'^ that 
J/(z,t)=E^ 



h-t 



dv e-"«^ A • P/(Z*'^ t + v)+ e-^^^'^-*) J/(4'!,, h) . 



(4.36) 
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With h = s, 

\Jf{z,t)-Jf{z,s)\<E^ 



e-^^^'X- Pf{Zl^\t + v)dv + \e 



^X(is~t) 



Jf{Zl%s)-Jf{z,s) 



s-t 



K{s-t) 



t,z 

s— 1> ' 



Jf{z,s) 



+ 



jf{z,s)y 

(4.37) 



In what follows we will bound the terms on the right-hand-side of this inequality. Since the 
condition (14.340 holds, Lemma l4. II applies. As a result it follows from ()4.1ip that 



{\Jf{Z';l„s) - Jf{z,s)\} <EE^ - .1} , 

Using the estimate in (j4.32p we have that 

y-1 



(4.38) 



EQ{Pf(Zi'M + v)} 



yi,{dy)E'^lfl^ + qtW 



y 



,t + v 



< I y„(dy)D(l + -E'i{\Zt'^\}+qt+J-^ 

'1+ \ y y 



< DU + E'^{\Z'/\} + {^ + l)qt+, 



< D{^ + —{C + l)+C{l + \z\ 



(4.39) 



To obtain the last inequality we use the estimation ()4.22p and the fact that qt+v < for 
V £ [0,s — t]. On the other hand, from (j4.33p . we have that 

\Jf{z,s)\<E{l + \z\). (4.40) 

In the inequalities above, the constants E, D and E are as in Lemma l4.ll and Remark 14. 2[ 
Now, using (|T38]1 . (|T39]) . (|O0]) and the inequalities 



e-^^'' < 1, and 1 - e'^^^^"*) <X^{s-t), 
we can bound (14.370 as follows: 
\Jfiz,t) - Jf{z,s)\ 



(4.41) 



< (^e + ^(e + l) + C(l + |z|)) ^s-t)+EE'i[\Zl'l,-z\] + X^E (1 + |z|) {s 



< DX [C + —i^ + l) + C{l + \z\) ] is-t)+E-C{l + \z\) is-t)2+XCE {l + \z\ 



< F {I + \z\) {s - t) 



-t) 

{s-t) 
(4.42) 



where -F is a positive constant only depending on A, 1^, T and Mf. To obtain the second 
inequality, we use the moment estimates (j4.23p . To obtain the last inequality, we use the fact 
that s-t<T. □ 
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In the following proposition we show that Jf satisfies a parabolic partial differential equa- 
tion. 

Proposition 4.1. Assume function / : M x [0,r] ^ satisfies the following condition 

\f{z,t)- f{z,s)\<D\z-z\+F{l + \z\)\s-t\"2, {z,t),{z,s)GRx [0,T], (4.43) 

IR+ is the unique classical 



in which D and F are constants, then the function Jf : M x [0, T] - 
solution, i.e. Jf G C^'^, of 

A{t)Jf{z, t) - XiJf{z, t) + X- Pf{z, t) + ^^Jf{z, t) = 
Jf{z,T) = {C-{z-K^))+. 



(4.44) 
(4.45) 



Proof. It is clear that Jf satisfies the terminal condition. For any point {z,t) G M x [0,T], let 
us take a rectangle R = [zi, Z2] x [0, T], so that {z, t) G R. Denote the parabolic boundary of R 
by doR := dR — [zi, Z2] x {0}. Consider the following parabolic partial differential equation 



d 

A{t)u(z, t) - X^uiz, t) + X- Pf(z, t) + —u{z, t) =0 

at 

u{z,t) = J f{z,t), on doR. 



(4.46) 
(4.47) 



Because of the condition (j4.43p . z — > f{z,t) is Lipschitz in its first variable uniformly in the 
second variable, it follows from Lemmas 14.11 and 14.31 that z — > Jf{z, t) is Lipschitz and t — > 
Jf{z,t) is Holder continuous. As a result J/(-, •) is a continuous function on M x M4.. 
On the other hand, for {z,t), (z, s) G R, it follows from the condition (j4.43p that 



\Pfiz,t)-Pf{z,s)\ < 



y y 



f {- + Qs- — -,s 

.y y 



yu{dy) 



< 



D\z - z\+ D\qt - qs\\y - l\ + F {y + \ z\ + qt\y - 1|) |s - t| 2 i^{dy) 



-rT 



< D\z - z\ + D{(, + 1) 



T 



e'^'du 



< D\z-z\+Fil + \z\)\s-t\2. 



+ F{i + qt{i + l) + \z\) \s-t\^ 

(4.48) 



in which F only depends on T and ^. Since i? is a bounded domain, the factor 1 + l^l in (j4.48p is 
bounded in i?, so 2: ^ Pf{z,t) is Lipschitz and t Pf{z,t) i s Holder, uniform ly with respect 
to the other variable. Now by Theorem 5.2 in Chapter 6 of Friedman ( 19751 ). the parabolic 
partial differential equation (j4.46p and (j4.47p has a unique classical solution in the bounded 
domain R. Moreover this solution can be represented by 



u{z, t) 



!^e-^^^JfiZt'^,r + t) + £ e-^«^A • P/(Z*'^ t + s) ds^ 



Jfiz,t), 



in which the exit time r = infgg[o,T-t] 

= zi or Z2} A (T — t). The second equality follows 
from the definition of the operator J in ()4.ip and the strong Markov property of Z*'^. 
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So far we have shown that Jf agrees with the unique classical solution of (j4.46p and (j4.47p 
in the bounded domain R. Since this statement holds for arbitrary R, it is clear that Jf is a 
solution of the parabolic partial differential equation (j4.44p and (j4.45p for all (z,t) gR xJO, T ]. 
The uniqueness of the solution follows from Corollary 4.4 in Chapter 6 in Friedman ( 19751 ). 
since the coefficients of the derivative operators in (j2.16p satisfy linear and quadratic growth 
conditions respectively. □ 

4.2 Properties of the Sequence of Functions Defined in ( 12.101) 

Our first goal is to prove z — > Vn{z,t) is Lipschitz and t — > Vn{z,t) is Holder continuous for all 
n. To this end we will apply Lemmas 14.11 and 14. 3[ To be able to apply the latter lemma we 
need to show that 

Mn= sup {?;„(0,t)} < oo, forn>0. (4.49) 

tG[0,T] 

In the next lemma, we will dominate the sequence of constants {Mn)n>o by a universal constant 
Moo, which depends only on T. 

Lemma 4.4. Let us define the sequence of constants {Mn)n>Q o.s in (f^.^g[ j, then 

A U a B 

Mn<Moo = + + i^i < oo, 

1 — a 1 — a ^ 

in which the constants U , B and a are defined in Lemma \4-3\ 
Proof. When n = 0, by the definition of fo(-, •) in (j2.10p . we have 

Mo = sup vo{0,t) = (C • (0 - i^i))+ < Ki, 

te[o,T] 

in which the last inequality is saturated when = —1. It follows from Lemma 14.21 that 

Mn+i <U + a (^Mn + ^) . forn > 0, (4.50) 
in which a < 1. It can be proven by induction that 

M„ < [/ ^^a* - a"^ +a ^^a^ - a"^ ^ + a"i^i, forn > 0. (4.51) 

Since U, B and ^ are positive constants and < a < 1, it is clear from ()4.5ip that 

U a B 

Mn < + + i^l = Moo < OO. 

1 — a 1 — a ^ 

□ 

Lemma 4.5. Let {vn{-, •))n>o be as in 112. 10\) . We have that 

\vniz,t) - Vniz,t)\ <\z - z\ , Z, 2 G M (4.52) 

and 

\Vn{z,t)-Vn{z,s)\<Fn{l + \z\){s-t)^, 0<t<S<T, (4.53) 

in which Fn are all finite constants depending on T. 
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Proof. ^From the definition of fo(-, •) in (j2.10p . we have 



\vo{z,t)-voCz,t)\ = \ {C-iz-K^))+ -{C-Cz-K^))+\ <\z-~z\. (4.54) 

Now, the inequahty (j4.52p follows from induction and Lemma l4.1i On the other hand, (j4.53p 
holds as a result of Lemma 14.31 which we can apply as a result of Lemma 14.41 □ 



As a corollary of Remark 14.21 and Lemma 14.41 we can show that t))„>o satisfies a 

linear growth condition in the z-variable, uniformly in the t-variable. This will be used to show 
that this sequence has a limit. 

Corollary 4.1. For any n> 0, 

Vn{z,t)<M^ + \z\=L{z), G M X [0,r]. (4.55) 

Proof. An induction argument using the inequality ()4.54p , Lemma 14.11 and Remark 14.21 gives 

Vniz,t) < Mn + \z\, for n > 0. 
Now, the result follows from Lemma |4.4[ □ 



As a result of Corollarv 14. H next we show that, for a fixed [z,t) G M x [0,T], the sequence 
{vn{z,t)}n>o is a Cauchy sequence. 

Lemma 4.6. For any (z, t) G M x [0, T] and n,m> 0. 



.i=0 



+ 2\z\B"', (4.56) 



where A = 1 — e ^^^'^ ^\ B = 1 — e ^^'^ and C is the same constant as in ^.2^ . 



Proof. We will prove the estimation (j4.56p by induction on m. When m = 0, it follows from 
Corollary lO that 

\vn{z,t) -voiz,t)\ < 2Moo + 2|z|. 



It is clear that (j4.56p is satisfied in this case. Assuming (j4.56p holds for m case, we will show 
that it holds when we replace m by m + 1. From the definition of {vn{-, ■)}n>o, we have 



\Vn+m+liz,t)-Vm+l{z,t)\ < E« <| ^ ds e~^^' X ■ \PVn+m{Z'/ , t + s) - PVm{Z'/ , t + s) 
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In the right hand side of above inequahty, the induction assumption gives us 

\PVn+m{Zl'',t + S)- PVm{Zl'',t + s)\ 



< 



' qt+s— — ,t + S ] - Vm 



y y 

< 2^A/oo ( 1 - e-^«(^-*-^)' 
+2^ 



s y ~ 1 

+ qt+s ,t + s 



yu{dy) 



1 ^ + 1 



+2 



t,z< 



+ qt+s — 
R+ V y y 



i=0 

y-1 ' 



1 - e 



-X{T-t-s) 



< 2^Moo(l - e 



-X£,{T-t)Yn 



+2C 



1 ^ + 1 



^(l_e-A(T-t)\ /^_^-A5(T-i) 



1 - e-^(^-*) 



j=0 



+2|Z*-|(l-e-^(^-))'" + -|(e + l) (l 



-A(T-t) 



(4.57) 



In (|4.57p . the third inequaUty follows, because qt+s < yp, and for m > 1 



m— 1 



m— 1 



-A(T-t-s) 



1 - e" 



-\^{T-t-s) 



El- 

1=0 

since s > 0. 

On the other hand, from (|4.7p . we have 



< 



i=0 



1 - e 



where E^Qll^,]} = E^{|Z*'°|} < C from ISM)- Therefore we have 

E^{\Zl''\} <\z\e^'' + C. (4.58) 
Taking expectation on both side of ()4.57p and plugging ()4.58p back into (I4.57h . we have 

\PVn+m{Zl'M + S)- PVm{Zi'\ t + s)\ 



< 2^MooA'" + 2^ 



1 e + i^c 



rT e e 



,i=0 



+2\z\ef''B"^ + 2 ( — (^ + 1) + C ) S™. 



(4.59) 



Multiplying both sides of (j4.59p with e •^^'^A and integrating with respect to s over [0,T — t], and 
using the identity // — = —A, we obtain the inequality (j4.56p with m replaced by m + 1. □ 

As a result of the previous lemma we can define the pointwise limit for the sequence 

(Vni-, ■))n>Q- 

Vociz, t) = lim Vn{z, t), {z, t) gRx [0, T]. (4.60) 

n>0 

Moreover, as a corollary of Lemma 14.61 we have 
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Corollary 4.2. For any compact domain D C M, Vn{z,t) converges uniformly to Voo{z,t) for 
{z,t) gV X [0,T]. Moreover, 

Mz,t)-vn(z,t)\ < Mr, (l - e-^^^^-'^^y , (4.61) 

where M-p is a constant depending on D and rj = max{^, 1}. 

Proof. Observing that the right hand side of (|4.56|) is independent of n and \z\ is uniformly 

bounded in D, the result follows from Lemma 14.61 □ 



In the following, we will begin to study properties of Vooi', •)• 
Lemma 4.7. The function Voo is a fixed point of the operator J. 
Proof. For any s G [0,T — t], 



\Zs\ \y-M 

Moo -\ h qt+s 

y y 



yu{dy) 



< eMoo + ^(e + l) + C(l + |z|). (4.62) 

As a result, we have 

Voo{z,t) = limvn+i{z,t) 
n>0 

= |e-^«(^-*) (C • [z'/_, - Ki)y + e-^«^A • (Plim ^;„)(Z*'^ t + s) ds^ 

= Jvoo{z,t). (4.63) 

The third equality follows by applying dominated convergence theorem three times. We can 
use the dominated convergence theorem due to Corollary 14.11 and ()4.62p . □ 

Using Lemmas 14.51 and ()4.6ip . we can show z — > Voo{z,t) is Lipschitz continuous and t — > 
Voo{z,t) is Holder continuous. 

Lemma 4.8. Voo{-,-) satisfies 

\voo{z,t)-Voo{z,t)\<\z-2\, for{z,t),{2,t)eRx[0,T]. (4.64) 

Proof. For fixed z and z, let us choose a compact domain V^^^ ^ so that z,z G Vz^z- Then 
we have 

\Voo{z,t) - Voo{z,t)\ < \Voo{z,t) - Vn{z,t)\ + \Vniz,t) - ■U„(z,t)| + \Vn{z,t) - Voc{z,t)\ 

< 2(1- e~^'i'^^'*'>Y Mr,^ . + \z-z\. (4.65) 



In order to obtain the last inequality, we use Lemmas 14.51 and Corollary 14.21 Since n in the 
second inequality in (I4.65P is arbitrary, the result follows. □ 
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Corollary 4.3. Voo{-,-) satisfies 

\voo{z,t)-Voc{z,s)\ <Foo(l + |z|)|t-s|5, (4.66) 
in which constant Foo < oo . 

Proof. This is a direct application of Lemmas 14.31 14.41 and 14.81 Note that Lemma 14.41 is needed 
to show that supjg^ 'p]{foo(0, t)} < cxo, which is required by Lemma [43l □ 

4.3 Proof of Theorem 

(i) This is a direct consequence of Lemma WM, which shows that the sequence {vn{z,t)}n>o is 
a Cauchy sequence. 

(ii) See Corollary [Ol 

(iii) Using the inequalities (j4.52p and (j4.53p in Lemma 14.51 we can apply Proposition 14. II to the 
function f = Vn- It indicates Jvn{-, ■) is the unique classical solution of the following equation 

A{t)JVn{z, t) - XCJVniz, t) + X ■ {PVn){z, t) + ^JVn{z, t) = (4.67) 

at 

Jvn{z,T) = iC ■ (z - K,))+ , 

for (z, t) G Mx [0, T]. By the definition of the sequence {vn{-, •))n>o in (|2.10p . we have Jvn{-, •) = 
Vn+ii', •)• So Vn+1 is the unique solution of (12.14p and (12.15p . 

(iv) Because of Lemma 14.81 and Corollary 14.31 we can apply Proposition 14.11 to the function 
/ = Voo- It shows Jvoo{-, •) is the unique classical solution of the following parabolic partial 
differential equation 

Ait)Jvoo{z, t) - X^Jvooiz, t) + X- {Pvoo){z, t) + -rrJvooiz, t) = (4.68) 

ot 

Jvoo{z,T) = {C-{z-Ki))+, (4.69) 

However, Jvoo = Voo by Lemma HTH Therefore, Voo{-, •) is the unique classical solution of the 
integro-partial differential equation p.l7p and ()2.18p . □ 

Acknowledgment We are grateful to the anonymous associate editor and the two referees 
for detailed comments that helped us improve our paper. 
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Table 1: The approximated price for continuously averaged European type Asian options for a 
double exponential jump model. 

r = 0.15, = 100, T = 1, p = 0.6 and ?7i = r/2 = 25. Monte Carlo method uses 10® 
simulations and 10"^ time steps. "C - P" is the difference between our approximated cah and 
put option prices. "Parity" is the difference predicted by the put-call parity. Run times are in 
seconds. 



European Asian call option prices for a double exponential jump diffusion model 









Iteration Algorithm 


Monte Carlo (Call Option) 


a 


K2 


A 


Call Option (C) 


Put Option (P) 


C - P 


Parity 














Value 


Time 


Value 


Time 


Value 


Std. Err. 


Time 




90 


1 


15.419 


1.0 


0.012 


1.0 


15.407 


15.398 


15.410 


0.006 


913 




3 


15.457 


1.5 


0.045 


1.5 


15.412 


15.442 


0.007 


976 


0.1 


100 


1 


7.170 


1.0 


0.376 


1.0 


6.794 


6.791 


7.170 


0.006 


919 


3 


7.456 


1.5 


0.656 


1.6 


6.800 


7.439 


0.007 


987 




110 


1 

3 


1.702 
2.220 


1.0 
1.5 


3.520 
4.040 


1.0 
1.6 


-1.818 
-1.820 


-1.817 


1.697 
2.207 


0.004 
0.004 


906 
981 




90 


1 


15.699 


1.0 


0.292 


1.0 


15.407 


15.398 


15.686 


0.012 


908 




3 


15.802 


1.5 


0.390 


1.6 


15.412 


15.806 


0.012 


983 


0.2 


100 


1 


8.540 


1.0 


1.745 


1.0 


6.795 


6.791 


8.540 


0.010 


935 


3 


8.790 


1.5 


1.994 


1.6 


6.796 


8.784 


0.010 


996 




110 


1 

3 


3.723 
4.045 


1.0 
1.6 


5.541 
5.864 


1.0 
1.6 


-1.818 
-1.819 


-1.817 


3.721 
4.038 


0.007 
0.007 


921 
983 



Table 2: The approximated price for continuously averaged European type Asian options for 
normal jump diffusion model. 

r = 0.15, ^0 = 100, T = 1, A = 1, /2 = -0.1 and a = 0.3. Monte Carlo method uses 10® 
simulations and 10'^ time steps. "C - P" is the difference between our approximated call and 
put option prices. "Parity" is the difference predicted by the put-call parity. Run times are in 
seconds. 



European Asian call option with normal jump 






Iteration Algorithm 


Monte Carlo (Call Option) 


a 


K2 


Call Option 


Put Option 


C - P 


Parity 












Value 


Time 


Value 


Time 


Value 


Std. Err. 


Time 




90 


16.997 


0.5 


1.601 


0.5 


15.396 


15.398 


16.991 


0.014 


913 


0.1 


100 


10.062 


0.5 


3.272 


0.5 


6.789 


6.791 


10.046 


0.013 


910 




110 


4.836 


0.5 


6.653 


0.5 


-1.817 


-1.816 


4.834 


0.011 


915 




90 


17.346 


0.5 


1.950 


0.5 


15.396 


15.398 


17.339 


0.017 


919 


0.2 


100 


10.959 


0.5 


4.170 


0.5 


6.789 


6.791 


10.968 


0.015 


917 




110 


6.303 


0.5 


8.120 


0.5 


-1.817 


-1.816 


6.310 


0.012 


913 
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The parameters for Asian options in the fohowing three tables are the same as the parameters 
used in the 7th row in Table 1, i.e. r = 0.15, = 100, Ki = 0, K2 = 90, T = I, a = 0.2, 
X = 1, p = 0.6 and 7?i = r/2 = 25. 

Table 3: The convergence of the option price with respect to the truncation length of the 
numerical integral. 

As we introduced in (|3.4p . the integral term in (|2.12|) is approximated by the trapezoidal rule 
on an interval [^^mmi •^max 

] with Xmin = xq < xi < ■ ■ ■ < XL = Xmax- In this table, fixing the 
discretization, we study the convergence with respect to the length of the truncation interval 
[xmin,Xmax]- We choosc Xmin = " -^/^2 and x^ax = ^/vi- In (j3.3p . if the distribution F is the 
double exponential, when N is large, the probability that the random variable X be outside 
the interval [—N/r]2, N/iji] is very small (for example, when = 15 the probability is less than 
10-6). 





Convergence 


with respect to truncation 




N 


Can Option (C) 


Time 


Put Option (P) 


Time 


(C-P) 


- Parity 


5 


15.5832 


0.500 


0.2858 


0.516 




-0.1002 


8 


15.6953 


0.765 


0.2916 


0.797 




0.0061 


10 


15.6994 


0.969 


0.2921 


1.000 




0.0097 


12 


15.6995 


1.141 


0.2921 


1.187 




0.0098 


15 


15.6995 


1.391 


0.2921 


1.500 




0.0098^ 



"Because we fix the discretization of the numerical integral, the difference between the calculated value and 
predicted value in the last column doesn't seem to converge to 0. But as Aa; 0, the difference will converge to 
as we will see in the next Table. 

Table 4: The convergence of the option price with respect to the grid size of the numerical 
integral. 

In this table, we fix the truncation of the numerical integral as Xmin = ~ 10/^/2 and Xmax = 
10/?/i, we will show the convergence with respect to the number of grids L in the discretization 
of numerical integral in (j3.4p . Since the density of double exponential distribution has a cusp 
at zero, we choose an unequaly spaced grid here. The closer x to zero is, the finer the grid is. 
While fixing the truncation interval [xmim^max], the larger L is the finer the grid is (see Table 
4 for the notation). 





Converj 


j;ence with respect to discretization 




L 


CaU Option (C) 


Time 


Put Option (P) 


Time 


(C-P) 


- Parity 


200 


15.7295 


0.422 


0.2926 


0.422 




0.0393 


300 


15.7103 


0.578 


0.2923 


0.609 




0.0204 


400 


15.7034 


0.766 


0.2924 


0.797 




0.0134 


500 


15.6944 


0.969 


0.2921 


1.000 




0.0097 


600 


15.6968 


1.141 


0.2920 


1.172 




0.0072 


700 


15.6954 


1.344 


0.2920 


1.360 




0.0058 


800 


15.6943 


1.516 


0.2920 


1.562 




0.0047 
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Table 5: The convergence of the option price with respect to the grid sizes used in the finite 
difference scheme. 

In this table we fix Xmin = — lO/r/2 and Xmax = 10/?7i, L = 1000 (See Table 4 for the notation). 
Moreover, we fix Zmin = z — 0.5 and Zmax = z + 0.5 with z = [l — e~^'^) /{rT) — e~^'^ K2/ Sq 
defined in (|2.20|) . We will show the convergence with respect to time and space grid sizes that 
are used in implementing the finite difference scheme. 



Convergence with respect to grid sizes 


Number of Time Steps 


Number of Space Steps 


Call Option Price 


Changes 


Time 


10 


40 


15.7093 


n.a. 


0.438 


25 


100 


15.6929 


0.0164 


1.890 


50 


200 


15.688 


0.0049 


7.500 


100 


400 


15.6864 


0.0016 


29.406 
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